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Abstract. - We investigate a modified Anderson model to study the local density of states 
(LDOS) of a molecular wire adsorbed on a metal. Using a self-consistent mean-field type 
approach we find an exponential decay of the LDOS along the molecule. A repulsive on-site 
interaction on the molecule suppresses the tunneling and decreases the characteristic decay 
length. 



Introduction. - In the last decade, a lot of progress has been made in the theoretical and 
experimental investigation of nanostructures like quantum dots and quantum wires Jl| . These 
structures are produced by lithography. Using this very flexible technology, a large variety of 
different structures with different geometries has been made. Simple dots, multiple dots, dot 
arrays and the corresponding leads have been realized both using metals and semiconductors. 
A disadvantage of lithographic methods is the presence of fluctuations in shape and size of 
these structures; they are unavoidable because of errors in writing the mask and because of 
the stochastic nature of the etching process. 

The idea to use molecules or supramolecular structures as quantum wires and quantum 
dots has been around for quite a while, see Ref. Q for an early suggestion of a molecular 
rectifier and Refs. || [| for recent reviews. Powerful chemical synthesis methods are available 
that allow the production of atomic and molecular clusters with linear dimensions of up to 
5 nm, i.e., they are approaching the lowest linear dimensions of nanostructures produced 
by lithography. It is a fascinating idea that such clusters which would have custom-made 
electronic properties could be included in nanostructures and replace, say, a quantum dot. 
Different copies of these systems would be identical in structure and properties (e.g., there is 
no difference in the composition of different Ceo molecules). First steps in the direction of 
this 'bottom-up approach' are under way in a number of labs around the world. Scanning 
probe techniques have been used to arrange atoms on a substrate in arbitrary ways || ||. 
A particularly impressive example is the 'quantum corral' j?J, a circularly shaped structure 
of Fe atoms on a Cu(lll) surface that may be considered a precursor of a device. Carbon 
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nanotubes have been contacted and proven to exhibit Coulomb-blockade behavior feL g| . A 
new kind of transistor whose central element is a Ceo molecule has been proposed [|l(J. The 
theory of conduction through molecular wires has been studied by a number of groups (see, 
e.g., Refs. ]ll|, |l2|, O , fL4| , [l5)). Also, the thermal conduction through molecular wires has 
been investigated ]16[]. Recently, the coupling of a molecule to a metal has been studied both 
experimentally and theoretically [|l7| |l8[ . 

A central problem facing the use of molecules as electronic devices is the question of how 
to contact them in a reproducible way. An interesting answer to this question was given in a 
recent experiment Q in which a custom-designed "lander" molecule Q was adsorbed on a 
(111) copper surface. The molecules can be thought of as consisting of an aromatic platform 
and four spacer "legs" . The spacer legs keep the main board at such a distance from the metal 
surface that the aromatic part is electronically decoupled from the surface. 

The molecules self-assembled perpendicularly to a step on the surface, thereby forming a 
contact between the metal and the aromatic board (the length of the legs was chosen to be 
comparable to the step height). A scanning-tunneling microscope (STM) was used to scan 
along the surface and the molecule, and it was found that the tunneling conductance decayed 
in an exponential way along the molecule. 

In this paper, we would like to study a (relatively) simple model for the molecule at the 
metal step. In contrast to the quantum-chemical calculations provided in |l^] that take into 
account the precise structure of the lander molecule, we propose to use a modified Anderson 
model to describe this situation. The 'impurity' of the new model is spatially extended, 
i.e., contains internal degrees of freedom not present in the standard Anderson model. The 
aim will be to calculate the local density of states of an extended structure (the molecule) 
using a Green's function method. The molecule will be approximated by a (finite) one- or 
two-dimensional tight-binding lattice with a Hubbard-like interaction. Our goal is to gain 
a qualitative understanding of the physical problem and its phenomenology (as opposed to 
quantum-chemical calculations for specific molecules) to be able to propose and analyze more 
complicated (e.g., multiply-connected) structures in the future. 

Model. - The system we have in mind is shown in Fig. |l|. Our model is defined by 
a Hamiltonian that consists of three parts, H = H c + Hd + Ht- The metal to which the 
molecular wire is attached is described by non-interacting electrons in a conduction band 

H c = 6kC L C kcr , (1) 

kcr 

where k and a =|, j denote electron wave vectors and spins, respectively. The single-electron 
energy ek is measured with respect to the Fermi energy (ep = 0). The molecular wire is 
assumed to be a one-dimensional lattice of L atomic sites and is described within the Hubbard 
model 

L L 

H d=Yl Mt^a + u J2 n omi • ( 2 ) 

i,j=l j = l 

The parameters that enter here are the on-site electron-electron interaction U and the matrix 
elements Uj = e^ij —t5ij±i : where is the on-site energy of the single-particle level at each 
site and t is the hopping element between nearest neighbors. The coupling of the wire to the 
metal is assumed to be small so that we can describe it by the tunneling amplitude Vk of the 
electron in the state (k, a) to the first site j = 1: 

H T = -J2(Vuct l7 d la + h.c.) . (3) 

kcr 
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Fig. 1. - Left part: schematic view of the system. We describe the molecule by a tight-binding 
model, i.e., a discrete lattice of atomic sites. The coupling of the molecule to a metal consisting 
of non-interacting electrons in a conduction band is assumed to be small and characterized by the 
tunneling amplitude Vk- Right part: sketch of the energy spectrum of the model, (a) unperturbed 
density of states N C (E) of the metal, (b) - (d) density of states p(E) at a typical molecular site, 
(b) Isolated molecule (Vk = 0) without electron-electron interaction (U = 0). (c) Molecule without 
interaction coupled to the metal. The coupling broadens the discrete molecular levels, (d) Molecule 
in the presence of strong electron-electron interaction (U ^> t) coupled to the metal. In the half- rilled 
case (one electron per site), the interaction will lead to the formation of an upper Hubbard band. 



The model defined by H is a generalization of the Anderson impurity model: for L = 1, 
the molecule reduces to a localized level, and (||) is the hybridization term of the Anderson 
impurity model. 

The qualitative nature of the spectrum of the system is shown in the right part of Fig. |[ 
We characterize the metal by a model density of states N C (E) to be discussed below and shown 
schematically in (a). If the molecule is not coupled to the metal, V = 0, its spectrum consists 
of discrete levels as shown in (b). Coupling the molecule to the metal broadens the spectral 
lines as shown in (c). If we assume a filling of one electron per site, the presence of strong 
on-site interactions leads to the formation of a Mott-Hubbard gap as shown in (d). 

We will use the imaginary-time path-integral formalism at finite temperature fc#T = 1/13 
(see, e.g., Ref. pip. In this formalism, the Euclidean action is given by 



f dr lj2 c LdrC ka + J2 d l d - d ^ + H ) ' 

J ° \ka ja J 



(4) 



and we want to calculate the thermal Green's function 

g(ja,ka'-r) = -(T T dU0)d ka ,(r) ) (5) 



and its Fourier transform defined by G(ja,k(j' ';iu n ) = Jq dr e l0Jnl 'Q(ja^ka' ';r), where u n = 
(2n + 1)tt/P with n integer are the Matsubara frequencies. 

First, we consider the non-interacting case, U = 0. In this case, the Euclidean action (|4|) is 
quadratic both in ^(cj^) and dj (T (dj (T ), and it follows that G(ja, ko' \ ico n ) = 5 aa 'G(j, k; iu n ), 

Q~ x ^,k\iuj n ) = g^d.k^iun) - E c (j, k\iu n ) , (6) 

or more explicitly 



(7) 
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The unperturbed (non-interacting, isolated) Green's function is given by 

Gq^ {j-> fc; iv n ) = iu n - tjk • 



(8) 



In Eqs. (||) and ([?]), the coupling to the metal manifests its effect on the wire through the 
self-energy 



F 

E c (j, k\iu n ) = SjiSki I 

J — c 



de r(e) 



2tt iu n — e 

The effect of the metal on the molecular wire is given by the function 

r(£) = -2]T|F k | 2 ImG?(k,£), 



(9) 



(10) 



where G^(k, E) is the retarded Green's function for the (unperturbed) metal. In the simplest 
case in which the energy-dependence of the tunneling amplitude Vk is not too large, Vk ~ V, 
( p7o| ) can be further reduced to an expression directly proportional to the density of states in 
the metal N C (E): T(E) = 2ttV 2 N c (E). 

The local density of states (LDOS) at the jth site is then given by the analytic continuation 
of the thermal Green's function. 



—lmQ(jJ]iu n 

7T 



E + i0 + ) 



(11) 



The factor 2 accounts for the contributions from the two spin components. 

We now turn to the opposite limit in which the electron-electron interaction on the molecule 
is very strong; U ^> t ]2^ ]. In this case, exact solutions are only available in special cases; 
the isolated molecule (Vk = 0), i.e, the usual Hubbard model Q and the single-site molecule 
[L = 1), i.e., the usual Anderson impurity model p4| . There are many approximation methods 
for the Hubbard model or the Anderson model (see Ref. |25| for a recent review), and we will 
adopt a self-consistent mean- field approximation |2l) . 

We start with the strong-Z7 limit of Hd in (||), i.e., the t—J Hamiltonian (2(J 



H d -^2Uj4 a d j(T - jj^Ujid^d^ - dl^tjkidj^ - d n d ki ) 



(12) 



ijk 



within the reduced Hilbert space without doubly occupied sites. The second term in ( |l^ ) 
can be made quadratic in dj C7 (d}- (7 ) by means of the Hubbard- Stratonovich transformation 
introducing the auxiliary field Aj to get the molecular part (corresponding to Hd) in the 
Euclidean action (0) 



qE _ qE 



Jo 



Sijd T + UjT3 



d 4 



(13) 



where rs is the Pauli matrix and Sj[ = Jq dr U^2- |Aj| 2 . In Eq. (|l3|), we have introduced a 
two-component spinor representation 



do 



4i 







A* 



(14) 



In this representation, the Green's function in Eq. (||) also has a matrix form ^(z,j;r) 



T r d i (0)dj(r) ). After integrating out the fields c kcr and c^, the Euclidean action can be 
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written as 



drdr' dt(r)S- 1 (^,i;r-r , )d,(r , ) 



where the Green's function Q is given by the Dyson equation 

<? _1 (j, k; iu> n ) = G^U, k; iu n ) - S c (j, k; iv n ) - E(j, fc; iw„) , 



(15) 



(16) 



with the new self-energy term E(z, j;r) = — 5(r)A^(r). Here (5o m @ and E c in Eq. (^|) have 
been extended to matrix form in a trivial way: Go = Goto and E c = E c fo, where ro is the 
identity matrix. 

So far no approximation has been made and the formal expression for the Euclidean action 
in Eq. (15) is exact. The interaction effect is correctly incorporated through the self-energy 
term E (whereas the effect of the metal is again manifested by E c ) as long as the integration 
over the field Aj is performed properly. At this point, we make our main approximation and 
neglect the fluctuations in the field Aj . We first assume a particular realization of the field Aj 
in the Dyson equation ( |l6| ) and then determine Aj within the stationary-phase approximation. 
Iterating this procedure allows us to determine Aj and the Green's function self-consistently. 
Moreover, the self-consistency equation 



A,; 



1 V- 



(17) 



corresponding to this procedure, gives us a physical interpretation of the field Aj: Aj measures 
the spin-singlet correlation between the jth site and its neighboring sites. It is clear that in 
the limit U — > oo, Aj w 2t/U unless j is too close to the boundaries. Finally, the LDOS 
can be obtained from the diagonal components of the retarded Green's function G R (j, k; E) = 
Q(j,k;iu n — > E + i0+): 



-— Im 

7T 



G? 2 (j,j; 



-E) 



(18) 



Results. - We will now report specific results on the LDOS of the molecular wire and 
discuss their physical implications. For the density of states (DOS) in the metal, we adopt the 
model functional form N C (E) = N C (0)(1 - E 2 /W?). The parameters iV c (0) and W c specify the 
DOS at the Fermi energy and the width of the conduction band, respectively. The self-energy 
contribution (||) from the coupling to the metal is then given by 



£ c (j, k; E) = S n S kl V 2 [a c {E) - iwN c (E)} 



with a c (E) defined by 



a c (E) = N c (0) 



2E 



1 



In 



E/W c + 1 



E/W c - 1 



(19) 



(20) 



All results shown in this work will be generated using the values €d = — 4t, V = O.lt, W c = 40t, 
and N c (0) — 10/t. The other parameters will be specified as needed. We also remark that we 
only investigate the lower band, i.e., the band below the Mott-Hubbard gap, in the interacting 
case. 

In Fig. we plot typical LDOS curves at three different sites (j = 1, 3, and 5) on a molecule 
with size L = 10 and on-site interaction U = lOt. The structures get sharper as U increases 
[ p7| because tunneling is blocked, see the discussion below. At a given site, the LDOS decays 
rapidly outside the band. 
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j=i 

j=3 
j=5 



-5 -4 
E/t 



Fig. 2. - Local density of states at site j for L = 10 and U = lOt. Solid line: 
(shifted upwards by 0.3), dashed line: j = 5 (shifted upwards by 0.6). 



1, dotted line: j — 3 



Figure || shows the spatial dependence of the LDOS along molecules with sizes L = 10 and 
20, respectively, both for the non- interacting (open circle) and interacting (solid circle) case. 
The LDOS has been calculated at the fixed energy E = — 2t, slightly above the upper edge of 
the band (remember that at the parameter values chosen below Eq. (20), E = — 2t is the upper 
edge of the band for an infinite tight-binding chain, i.e., Hd with L = oo and U = in Eq. @). 
The LDOS decays exponentially with distance from the metal- molecule contact, except for the 
region close to the boundaries. This reproduces the exponential decay of the STM tunneling 
conductance in Ref. ]l9[ . It also agrees with the expectation that the electron wave function 
induced by the metal to the molecule should decay in an exponential way. Another interesting 
conclusion that can be drawn from Fig. ^ is the influence of the interaction on the suppression 
of the LDOS. The characteristic decay length decreases with the interaction U. This can be 
interpreted by saying that the strong repulsion by the electrons sitting on (singly-occupied) 
sites tends to block tunneling events from the metal to the molecule. This effectively leads to 
the suppression of the "hybridization" term Ht, Eq. ([|). The suppression of the LDOS by 
interaction effects does not depend on the choice E = — 2£, however, the spatial scale would be 
different for different positions of the molecular level. The influence of repulsive interactions 
was also studied in (2^] , where it is shown that the suppression is lower for repulsive interactions 
than for a dimerization leading to the same value of the gap. 

In conclusion, we have calculated the local density of states of a molecule adsorbed to 
a metal surface using a new Anderson-type model. The LDOS decays along the molecule 
in an exponential way (like in the experiment ]l9|) and is suppressed by a local interaction 
on the molecule. Our formalism provides a qualitative understanding of this and similar 
experiments. It is able to include local interaction effects and can be used to treat more 
complicated geometries that will be investigated in future experiments. 



We would like to acknowledge discussions with T. Jung and A. Baratoff and thank 
J. Gimzewski for providing a preprint of Ref. [[19|. 
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